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Abstract 

The two major functions in galaxy evolution that we would like to measure are the 
stellar populations in galaxies, and their time-derivative, the star formation rate. 
Especially at redshifts above 1, both of these measures are benefiting greatly from 
wide-area infrared observations. Several space-based and ground-based examples will 
be discussed. 

1 Understanding Galaxy Evolution with Deep Fields 

Since this is a conference talk rather than a journal paper, I'll start on a philosophical note. Let's reflect 
on some of the good fortune observational cosmologists enjoy. The task of understanding the contents 
of the universe, (let alone their evolution), is so audacious, that it is remarkable we can even start doing 
the job. One reason we can is that gravity helpfully imposes some clear organization on how matter 
is arranged. Luminous mass appears (at most wavelengths) to be strikingly organized into stars and 
galaxies. Biologists rely heavily on their crucial organizing concepts of cell, organism and species. As 
astronomers, we see the world clearly arranged into self-gravitating planets, stars and galaxies. Although 
these categories span enormous ranges, for example in luminosity, they obey regularities and patterns. 
This is why we can make clear distinctions between what is or is not a star, and what is or is not a 
galaxy. We can, for example, distinguish individuals in these classes from bound clusters of them. These 
are two of the surprisingly easy categories which define our field. When we cannot use them, our job gets 
harder. Operationally, the star and galaxy concepts are closely linked, since we rarely think about one 
being present without the other. It is in fact very difficult to learn much about "galaxies" before they 
began forming many stars. I'll discuss one possible exception to that generalization at the end. 

The last half of the 20th century saw dramatic progress in understanding the evolution and even the 
formation of stars. The first half of the 21st will herald similar advances in understanding the evolution 
and formation of galaxies. These systems are sufficiently more complicated and distant than stars that 
empirical, observation- driven data must play a relatively larger role than theory. (Stellar structure 
and evolution are rarely influenced by external circumstances. Isolated stars, for example, hardly ever 
experience close encounters. Galaxies often do, and their motions are often dominated by dark matter 
we know little about.) 

Since galaxies have a wide range of (evolving) properties, we need to study large samples of them. 
Fortunately, they are so numerous that deep field surveys with large array detectors can provide such 
samples, often spanning a wide range of redshift simultaneously. Two of the most fundamental observables 
to derive from galaxy surveys are 1) their luminosity function (a measure of their population of formed 
stars) and 2) their rates of recent and ongoing star formation. In deep fields we can study how both of 
these evolve through cosmic time. 

Of the results obtained from the Hubble Deep Field (HDF), the most far-reaching is that the star 
formation rate is high at z ~ 3 ([p0[,|pl|), with a possible decline at z ~ 4. When combined with low- 

redshift data (e.g. Esl), there is evidence for a peak in the comoving SFR at z w 1 2, around a third 

of the Hubble Time. The apparent resemblance between this star formation rate curve and the evolution 
of quasars has been noted. 

For both of the main goals in the study of galaxy evolution-measuring the stellar population, and 
its first time derivative-observations with large-format (e. (/.million-pixel) infrared detectors are essential. 
I will illustrate this by reviewing several applications of infrared observations. Some of these were also 
discussed at a Tokyo conference of the Birth and Evolution of the Universe (]24|). Even though only two 
years have passed, so many advances have been made, that an update is due. 
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2 Detecting Evolving Stellar Populations Across Cosmic Time 



Deep multicolor imaging is our most powerful tool for understanding galaxy evolution. Broadband filter 
measurements can go very deep; the resulting spectral energy distributions provide estimates of the 
photometric redshifts of distant galaxies. The Hubble Space Telescope, with its superb spatial resolution, 
provides deep detections, along with detailed measurements of galaxy morphologies, in small fields. Large 
CCD cameras on ground-based telescopes provide complementary coverage of much larger areas (e.g. in 
PH , and the Sloan Digital Sky Survey.) The next step, intensive spectroscopic follow-up (e.g., has 
been rewarded with the detection of many hundreds of galaxies at redshifts of around 1 and 3. Most 
galaxies at z > 3 have been discovered using the powerful Lyman break/U-band drop out selection 
method (e.g., [0). This method, which is currently confined to z ~ 3 - 4, does not of itself provide an 
accurate confirmed redshift. Follow-up optical spectroscopy is still needed. The Lyman break method 
falters at z < 2.7 because at that redshift the Lyman limit has not yet cut out much of the continuum in 
the U-filter. Lyman break galaxy (LBG) surveys are not possible in the redshift 2 range. 

The intermediate redshift range around z ~ 2 appears to be the epoch during which the familiar 
components of modern galaxies assembled. A significant percentage of the universe's evolution (« 3 
Gyr), or 25% of its total age, occurred from z=l to 3, as opposed to the more distant and better-studied 
z=3-4 range, which covers « 0.5 Gyr, or only 5% of cosmic time. This intermediate epoch can be dubbed 
the "Bright Ages" , since it appears to be the time when most stars in the Universe formed (and most 
heavy elements were produced). 

However, our observational knowledge of galaxies in the crucial range of z = 1.5 — 2.5 is surprisingly 
sketchy. Due to severe limitations of optical search and confirmation methods, we are in the odd position 
of knowing more about the most distant galaxies than those closer. Continuum observations at A > 4000A 
are needed to assess the build-up of evolved stars. Similarly, the most useful emission lines for measuring 
gas photoionized by recently formed stars are in the rest-frame optical. Thus the most reliable measures 
of total stellar mass as well as current star formation rates (the first derivative of the former) are obtained 
in the near-infrared for redshifts in the Bright Ages. 



2.1 Near-Infrared Searches for Galaxies in the Bright Ages 

Some of the early successes in identifying non-active, non-lensed galaxies at high redshifts came from 
searches for strong emission lines with narrow-band filters. This method provides accurate redshifts, but 
selects from a limited range of redshifts and has yielded only a modest number of detections (e.g., p7[ , 
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Optical multicolor photometry alone works well at measuring the Balmer break up to z ~ 1. Infrared 
photometry makes it possible for us to identify Balmer break galaxies at 1.5 < z < 2.5. New large- format 
detectors are able to make sensitive surveys in the near-IR wavebands, especially with the key 2.2 /im 
(K) band. This allows us to measure the second strongest spectral feature in galaxies-the Balmer break, 
from around 4000A to 3650A. This blue-side drop is nearly always strong in populations of young or old 
stars (although its wavelength shifts slightly redward in older galaxies). 

Figure 1 shows the V, I, J and K bandpasses superposed on the model spectrum of a galaxy at z=1.5. 
To reach z=2.5 requires good photometry in the K filter, so that at least one point in the spectral energy 
distribution is cleanly on the red side of the Balmer break. Two-color plots are useful for identifying 
galaxies likely to lie in the Bright Ages epoch. An example is shown in Figure 2. Models with a variety 
of star-formation histories are shown for galaxies with < z < 1 by the crosses; models with 1 < z < 2 
by the diamonds, and 2 < z < 3 with the small dots. In the astro-ph version I color-coded these model 
predictions based on their V-I colors. The open squares show galaxy photometry from one of our deep 
fields imaged on the Lick 3-m telescope. 

The 3600/4000 A break is less strong than the Lyman break. Nonetheless, by fitting the full galaxy 
multiwavelength spectral energy distribution to redshifted galaxy templates we can estimate its photo- 
metric redshift to an accuracy of one or two tenths in z (e.g., |17 , 0, jjj, |ll|). Infrared flux points 



greatly improve the accuracy of the method To demonstrate the power of this method, we present 
a six band (B through K) picture of an very red object |^|. Figure 4 shows the photometric redshift 
fit to its BVIJHK spectral energy distribution, using Hyperz (Q). The infrared bands were crucial not 
only in identifying this as an interesting object, but allowing a fit to data beyond the 4000A break. The 
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Figure 1: The VIJK filter set superposed on typical galaxy spectra at z=1.5. The representative Spectral 
Energy Distribution for the galaxy is from the GISSEL spectral synthesis models. The Balmer break 
produces a very red I- J color. As the redshift increases above 1.5, the shift in the Balmer break makes 
J-K redder than 2. 




Figure 2: The the I-J,J-K plane for model galaxies at redshifts from to 5. In the astro-ph version, 
additional V-I information is encoded in color (with the bluest values shown as the bluest points). The 
colored points outline tracks for Bruzual and Chariot stellar evolution models of galaxies with initial 
formation epochs at high redshift which then have superposed more recent starbursts of varying strengths. 
The limited grid of starburst epochs (z=2.6, 2.8. 3.0, 3.5, 5) produces the several different locations of 
the "peeling back' in the two-color diagram. Note that the rest-frame UV light (shifted into the B+V 
bands) is extremely sensitive to the temporary presence of young stars, while the rest-frame optical (in 
the J and K bands) is insensitive to them. The open squares show our measurements of galaxies in one 
of the deep fields we observed with the Gemini 2-channel infrared camera on the Lick 3-m telescope. 
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Figure 3: Mosaic of reduced multicolor images of a 4 x 4 " box around a galaxy in the 53W002 field with 
photometric redshift z = 1.12. Wavelength increases from B at the left to K on the right; the object 
can qualify as an "Extremely Red Object" (ERO), with V — K > 6. The smooth round bulge becomes 
progressively more dominant at longer wavelengths, while at the shorter wavelengths the galaxy appears 
patchy with irregular arms. 
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Figure 4: BVIJK photometry we measured for an "extremely red" galaxy in the MTM0953+549 field. 
Note the sudden jump at the Balmer break, which makes a very large V-J color, even though the J-K 
is not unusually red. The best-fitting spectral template, at z=1.23 is overplotted in the top panel. The 
lower panel shows the probability distribution of acceptable redshifts fitted by the photometric redshift 
engine. The best fit implies the galaxy has My = —22.9. 



fit estimates its redshift, z=1.23 +- 0.25, and also that it is a 0.5 Gyr old burst of star formation with 
Ay = 1.2 mag of dust extinction. The infrared photometry provides a fairly good estimate of the total 
stellar mass. 

A lot of what we think we know about high-redshift galaxies rests on the limited statistics of photomet- 
ric redshifts in the Hubble Deep Fields, mostly HDF-N. The HDF-N has 150 galaxies with photometric 
redshifts in the Bright Ages (1.8 < z < 2.6). However, only 8 of these are spectroscopically confirmed 
after extensive efforts, and no Bright Ages galaxy candidates in HDF-S have confirming spectra yet. 
Major conclusions have been drawn from analysis of one 2.6 x 2.6 image of the Universe. This is par- 
ticularly dangerous since LBG surveys have found large field-to-field variations-up to a factor of four in 
surface density-due to large-scale structures, which also appear as spikes in redshift histograms. Since 
our view of galaxies depends strongly on which sight-line we observe, it is imperative to investigate galaxy 
evolution in several regions with similarly deep data. Many efforts, both ground- and space-based, are 
currently underway to do this, for example at Subaru ([§2|]), VLT (VIRMOS/DEEP) and NOAO (Deep 
Wide-Field Survey). 



4 



Figure 5: Our WFPC-2 continuum-subtracted Ha image of "Early- Type" Galaxy NGC 2328 from Glass- 
man and Malkan flip] , shown as a greyscale plot. The optical continuum is shown with overplotted 
contours. The field of view is roughly 6 x 10". 



3 Evolution of Star Formation is also well traced in the infrared 

A region's "current" star formation rate (SFR) is approximately proportional to its population of upper 
main sequence stars, because of their cosmically brief lifetimes. The presence of massive stars is manifested 
in several observables which are not strongly produced by old stellar populations: a) direct photospheric 
radiation, predominantly in the UV. At lower frequencies it extends into the optical, where it produces 
blue colors in a stellar population. The high-frequency tail extends above 1 Ryd, and includes a substantial 
fraction of ionizing photons, in the hottest stars. Much of this radiation is absorbed by gas and re-emitted 
as b) nebular emission lines, mostly in the optical and infrared. The most massive stars hardly live long 
enough to completely leave the dusty environments in which they were born. Thus, much of their 
remaining continuum emission may be absorbed by interstellar dust grains, which are warmed to c) re- 
emit the power thermally in the infrared. Finally, massive stars are the most violent sources of mechanical 
energy in the ISM, driving winds and supernovae. A fraction of this power accelerates electrons which 
then emit d) synchrotron radiation. 

After a great deal of work on the alternatives, it is Ha line emission which remains the "gold standard" 
of SFR measurements. It is equivalent to counting the ionizing continuum photons from the hottest 
(youngest) main-sequence stars at a wavelength that is relatively insensitive to absorption. The observable 
line luminosity is supposed to equal a constant times the aggregate SFR: L(Ha) = 1.1 X 10 41 erg/sec per 
Ms/year (T(J. 

By contrast, the other leading SFR indicator at high redshift-the ultraviolet continuum-is far more 
sensitive to dust reddening. The integrated star-forming luminosity of those LBGs detected needs to be 
corrected upward for extinction by a factor of 2 - 7 (| 32 , (TJ], and J38|). Furthermore, these 

LBGs were selected to be the least dusty ones. 

We can obtain a better idea of why the line and continuum SFR estimates disagree (methods a and b 
above) by looking closely (with HST) at nearby star-forming galaxies. The contours in Figure 5 show our 
HST imaging of the (continuum-subtracted) Ha emission in the small star- forming galaxy NGC 2328. 
The greyscale shows the optical continuum, which has a considerably different spatial distribution. Note 
the numerous contour peaks with little or no associated optical light. These HII region complexes with 
significant reddening account for the majority of the current star formation in this galaxy. Another even 
more dramatic example of a dust-enshrouded starburst is seen in the nearby galaxy NGC 5253 jyj. 

Although Ha gives a more complete census of recent star formation, it cannot be measured with CCD 
spectrographs for z > 0.5. Even the weaker and less reliable H/3 line is shifted out of the optical range at 
z > 1. Fortunately, they are still measurable in the near-infrared. Two methods-one ground-based and 
one space-based-have been successful. 



3.1 Narrow-Band Infrared Imaging of Emission-Line Galaxies 

Large-format multislit spectrographs are not yet available in the ncar-IR. This will change when instru- 
ments such as the NIRMOS (Near-infrared Multi-Object Spectrograph) begin working on the VLT, and 
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Figure 6: Extracted spectrum of a star-forming galaxy discovered in a deep NICMOS-grism parallel 
observation (Malkan et al. J23|). At a redshift of 1.49, this low-resolution spectrum shows both Ha and 
H/3 , as well as a broad blend of the two [OIII] doublet lines, the redshifted 5007 and 4959A emission. 



Flamingos (in multi-slit mode) on Gemini. For now, ground-based searches are restricted to either narrow 
spatial or spectral windows. The latter case is obtained with narrow-band interference filter imaging. To 
speed up the detections, there is often a pre-selection of special fields. These are supposed to have excess 
galaxies at the targeted redshift, because some object in the field is already known to have that z. Several 
dozen galaxies at redshifts above 2 have been discovered by detecting their Ha line emission (p7y, p8| , 
||,|Lq], p9[). More discoveries will be coming soon with the availability of IK x IK near-infrared array 
detectors. 



3.2 Grism Spectroscopy 

Slitless spectroscopy receives the full sky background from the entire spectral window at each point on 
the detector. It is therefore impractical with ground-based infrared telescopes. However, it has been 
demonstrated to work very well in space, where the backgrounds are much lower. In particular, the 
NICMOS Camera 3 on HST has an objective grism that disperses the spectra of all objects within its 
50" field-of-view. Since it was able to operate simultaneously with other HST instruments, it obtained 
a substantial number of deep parallel pointings on random fields. The largest set of these were reduced 
and analyzed by McCarthy et al. |3(|. The exposures of a few orbits had a 1 in 3 chance of detecting a 
line-emitting object in the field; in the longer parallel observations, these odds increase to nearly 100%. In 
addition, a single very deep NICMOS grism field was observed during a 3-day pointing in the Continuous 
Viewing Zone (|2|). 

The standard assumption in deep near-IR grism surveys of random fields is that any single emission 
line detected is redshifted Ha . In only one of the dozen cases tested so far has this assumption been 
shown to be incorrect; it is likely true better than 90% of the time. This allowed Yan et al. jl2| to 
estimate the global luminosity function of recent star formation in the 0.8 < z < 1.7 universe. 

In some cases, the difficult confirmation with optical spectroscopy is not even necessary because a 
second emission line is also detectable simultaneously in the NICMOS grism image. Figures 6 and 7 show 
deep grism spectra (observed wavelengths) in which emission lines other than Ha are detected. In Figure 

6 the additional lines on the blue side are Hp, and the [OIII] 5007/4959 doublet, at z e = 1.49. In Figure 

7 the only emission-line spectrum which does not show Ha is shown. This foreground (z e = 0.35) galaxy 
shows the strong emission lines of He I 1.083/J>7 1.093 ^m, [Fell] 1.257 and P/31.281 ^m ||. 

In its first incarnation, this camera was defocused by the distortion of the NICMOS dewar. This 
degraded the spatial resolution to about 0.4". The resulting spectral resolution was only R ~ 50. 
Hopkins et al. [Q used NICMOS Camera 3 for pointed observations of suspected high-redshift galaxies 
during a campaign when it was in focus. Both their observations and the parallel ones indicated that 
SFR rates estimated from Ha emission in galaxies at 0.8 < z < 1.7 are several times higher than those 
estimated from their UV continuum. 

These large numbers of strong Ha -emitting galaxies already indicate that search methods using 
shorter wavlengths are incomplete. Similar incompleteness in surveys of star formation rates has also 
been found in the local universe by Sakai et al. |35| . If short- wavelength surveys really do miss a substantial 
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Figure 7: Another star-forming galaxy in the same NICMOS parallel field |23|. This is the only grism 
spectrum in which none of the detected emission lines is Ha . The lower redshift, z e = 0.35, allows the 
detection of the He I I.O83/P7 1.093 /im blend and the [Fell] 1.257 and Pj3 1.281 jxm emission lines. 



Comparison of SFR Estimators in 
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Figure 8: Comparison of star formation rates estimated from [OII]3727 and Ha line emission, in z > 1 
galaxies discovered by the NICMOS grism. Typical uncertainties on the points are 20-30%. The open 
square indicates the line luminosities predicted by the local relations for an SFR of lOM^/year. All of 
these 9 galaxies observed at Keck by Hicks et al. [[12) have relatively weaker [Oil] emission than would be 
predicted by their Ha line, based on local spiral galaxies (which are described by the solid line). Instead, 
their [Oil] /Ha ratios are suggestive of additional internal dust absorption, which might be as large as 
Ay — 3 mag (shown by the dashed line). 



7 



in 



I 



\ 4 - _ 



O 



A 



O 



0.0 0.5 1.0 1.5 2.0 
(R-K) 

Figure 9: The ratio of star formation rates estimated by [OIII]5007 and UV continuum emission from 
10 Lyman break galaxies [Q. [OIII] for the open circles and asterisks was measured from narrow- band 
imaging; for the remaining points it comes from slit spectroscopy. The [OIII] line reveals several (up to 
seven) times more current star formation than does the UV continuum, especially in the redder galaxies. 
This is probably largely attributable to internal dust reddening. 



amount of obscured star formation activity, we should examine galaxies that are found in the near-IR, 
to see if they are indeed dustier than those discovered by optical surveys. In Figure 8 we compare the 
two emission-line indicators of current star formation rates-Ha and [OII]3727 luminosities. The nine 
data points represent the NICMOS-selected galaxies, with followup Keck optical spectra from Hicks 
et al. (l2). The solid line shows the observed average ratio in local spirals, which does not fit the 
NICMOS galaxies. In these IR-selected galaxies the [Oil] line systematically underestimates the SFR. 
This discrepancy between the two indicators cannot be predicted from the broadband colors, which are 
not often extremely red. Many of these galaxies must contain additional internal reddening around their 
star-forming regions. The dashed line shows an indicative amount, of Ay = 3 mag. If the ten galaxies we 
measured are representative of Ha selection in general, it produces samples of galaxies which a) would 
mostly be missed by optical search methods; and b) still need a significant upward correction for dust 
extinction. We need Ha and short-wavelength SFR measurements for more of the same galaxies to see 
how much overlap there is between those selected with optical and IR searches. 

3.3 Nebular Diagnostic Line Emission at Higher Redshift 

Ha is redshifted beyond the prime sensitivity of infrared detectors at z > 2.6. By unlucky coincidence, 
the rising thermal background makes it impossible to compare star formation rates from Ha with those 
estimated from the ultraviolet continuum in the LBGs. Thus the only way to make a direct comparison 
of SFR indicators (i.e., on the same galaxies) is by replacing Ha with another emission line that is bluer 
in the optical rest-frame, and less reliable. The advantage of using H/3 as this "silver standard" SFR 
estimator is that its intrinsic strength relative to Ha is known accurately from recombination theory J33| . 
The [OIII] 5007 and 4959A doublet, however, is usually two or three times stronger, and is not so affected 
by underlying stellar absorption. In Figure 9 we compare the star formation rates predicted by [OIII] 
and UV continuum emission in eleven LBG's at z ~ 3. A few galaxies lie near the bottom of the graph, 
consistent with 1-to-l agreement. However in most of them, the [OIII] line emission predicts a higher 
SFR by a factor of 3 to 7 times. This is again consistent with the amounts of internal dust reddening 
that have been inferred from other observations. 
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Figure 10: Model predictions of the Diffuse Infrared Background Emission from galaxies, compared with 
observations, from Malkan and Stecker [^6|. The "best guess" baseline evolutionary model shown by the 
middle line, assumes that infrared luminosities of galaxies increase as (l + z) 3 from z=0 back to z=2. The 
most rapid evolution model, shown by the upper ("Fast") line, assumes luminosity evolution as (1 + z) A 
back to z=1.4, and is also marginally consistent with all the observations. The lower limits are set by 
deep number counts obtained with the Infrared Space Observatory. The two lines with upper limits are 
set by the absence of high-energy absorption in the gamma-ray spectra of distant blazars, from Stecker 
and DeJong. Most of these strong observational constraints became available after the MS models were 
first published. 



4 The Promise of Space Imaging Surveys in the Thermal In- 
frared 

The thermal dust emission, (method "c" above) can estimate the high-redshift SFR with hardly any 
sensitivity to extinction. Instead of measuring the primary signatures of hot young stars- their UV 
continuum and resulting ionized gas-we can measure this emission after it is reprocessed by dust, since 
that is where much (probably most) of the energy ends up before escaping the galaxy. This secondary 
thermal continuum is emitted mostly at 25-60 fim. Since this continuum emission can dominate the 
entire bolometric luminosity of a dusty star-forming galaxy, it is now detectable out to high redshifts. 

First, it can be detected globally, by wide-area measurements of the Diffuse Infrared Background 
(DIRB) radiation. New observations continue to define the DIRB more and more accurately, over a wide 
wavelength range (from 2-500 ^m). It is surprisingly bright, and requires strong cosmic evolution in the 
population of dusty galaxies at z ~ 1 (e.g., p5||). "Backwards evolution" models, which are based on our 
empirical knowledge of infrared emission from local galaxies, have also improved. Figure 10 illustrates 
a recent comparison between these models and DIRB observations, from Malkan and Stecker J26|. Two 
plausible evolutionary scenarios are shown - a baseline model with pure luminosity evolution going as 
(1 + z) 3 , and a "Fast" evolution scenario that reaches the same maximum luminosities at an earlier z, 1.4 
rather than 2.0. As shown in Figure 11, this latter scenario is somewhat more consistent with independent 
estimates of the cosmic SFR evolution, compiled by Rowan- Robinson ||34) . 

Upcoming deep IR surveys from space will soon resolve this background. The key is surveying enough 
area with sufficient sensitivity to detect "typical" galaxies at the desired redshifts. "Typical" can be 
defined by the characteristic knee in the galaxy luminosity function, which is well defined at z=0, and 
appears to be more luminous at high redshifts |4C|] . Although such "Z*" galaxies are not very numerous 
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Figure 11: Comparison of the luminosity evolution model from Malkan and Stecker (1998 and 2001, 
p5[ , p6[) with estimates of the global evolution of star formation rates. The baseline evolution, and 
rapid evolution models of MS are shown by the light and dark solid lines, respectively. Both assume no 
evolution from z of 2 up to 5. Error bars are estimates from Rowan- Robinson p4| ] based on a variety of 
observational methods, with substantial corrections for dust extinction applied. 



in a standard Schechter luminosity function, galaxies of L* or brighter do produce more than 16% of the 
total integrated galactic luminosity. And in flux-limited samples, where the number of galaxies detected 
gets weighted like L 3 / 2 , half of the total luminosity will come from galaxies of 0.7 £* and brighter. Figure 
12 shows the Malkan and Stecker predictions for the observed fluxes of galaxies back to redshift 5, 
for 3 different wavelengths. Note that they do not fade much from z—0.5 to 2 because of their strong 
luminosity evolution. 

SIRTF, in its "First Look Survey", and Astro-E (IRIS), in a wider sky survey, will have the sensitivity 
needed to detect the typical £» galaxies at z > 1 which are expected to dominate the DIRB. As shown in 
Figure 13, the longest cosmological reach is provided by observations at 6 /im (with SIRTF) and 200 /im 
(with FIRST), since these wavelengths are redward of the two peaks in the spectra of galaxies. This 
results in "negative K-corrections" , the reverse of the usual situation in which higher redshifts shift a 
relatively fainter part of the galactic spectrum into our fixed observing bandpass. 



5 Damped Lya Absorption from Pregalactic Objects? 

Implicit in the previous discussion is the assumption that the "birth of a galaxy" is operationally equiv- 
alent to the epoch of its first major burst of star formation. This is usually the first chance we have 
to detect objects which will eventually evolve into modern galaxies. However, damped Lya absorption 
systems (DLAs), detected as neutral hydrogen troughs (Njji ~ 10 20_23 cm^ 2 ) in the spectra of high- 
z quasars, detect baryonic concentrations independent of star formation or nonstellar nuclear activity. 
DLAs contain most of the H I in the universe. 

Imaging searches for galaxies associated with DLAs are not new, but they have been almost impossible 
to carry out from the ground [19). They are best performed in the near-IR, where the DLA starlight 
should be relatively bright compared to the background quasar. Unfortunately, a large enough H I 
column density to be detected as a DLA is only expected to occur at small impact parameters through 
a galaxy (or protogalaxy) - within 10 or 20 kpc of its center (^J). This corresponds to a separation 
between the galaxy and the line-of-sight to the background quasar of less than 1 to 2". Detecting a faint 
galaxy this close to a bright source is exceedingly difficult from the ground. Adaptive optics is not yet 
the solution (beset by poor Strehl ratios). Proposed interferrometric techniques are more suitable for 
detecting a compact source at a known position angle with respect to the quasar (and still the required 
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Figure 12: Predicted brightness at infrared wavelengths of an galaxy, from Malkan and Stecker p6j, 
at 25, 60 and 200 /im. The dotted lines show the predicted brightnesses for galaxies 100 times above and 
below the characteristic L* knee of the luminosity function. There is a kink at z=2, beyond which the 
galaxy luminosity evolution is assumed to level off. 
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Figure 13: Predicted logarithmic number counts per unit redshift for galaxy luminosity evolution models 
in Malkan and Stecker p6|. The higher line, shown by dashes, illustrates the "Fast evolution" model. 
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Figure 14: Upper limits to the 1.6 /im continuum emission from starlight associated with Damped Lya 
Absorbers, from Colbert and Malkan (|J). Since they are all limits (except for the two plausible detections 
shown by asterisks), they are pointing upwards (towards less luminous galaxies). These upper limits would 
be a factor of two tighter for any galaxies more than 1" away from the quasar (ie., they understate the 
true sensitivity of the observations). The dashed, dotted and solid lines illustrate typical 5 — a upper 
limits obtained in our Cycle 7 program, and future possible NICMOS observations with somewhat longer, 
and much longer integration times. For twenty observed DLA systems, our HST coronographic imaging 
rules out the possibility of a galaxy having typically a luminosity of 0.7-1.0 L. t . 

dynamic range may be unattainable). Ground-based imaging can discover DLA galaxies only at large 
impact parameters, greater than 25 kpc (for H Q = 75). To do better, an extremely sharp and stable PSF 
is required, preferably near the diffraction limit. This means going into space. 

We therefore used the NICMOS-2 camera on the HST in coronographic imaging mode to search for 
near-IR continuum emission from starlight associated with two dozen DLAs. Our results are shown 
in Figure 14. Despite the approval of 95 targets, only two dozen quasar observations were attempted, 
and many were not as sensitive as they could have been. This low success rate was because: 1) the 
coronograph was not fully operational during the first half of Cycle 7, 2) SNAPSHOT proposals utilizing 
the coronograph were not scheduled for a few months after the coronograph became operational, and 
3) failure of the coronograph to autonomously acquire several of the quasars when they were scheduled 
resulted in a loss of those observations. Despite these difficulties, we detected two likely DLA galaxies of 
~ 0.7L*, and obtained upper limits for nearly twenty DLAs (j^])- Our 5<7 upper limits for the continuum 
fluxes of the DLAs are plotted in terms of versus redshift in Figure 14. Note that the vertical axis has 
logarithmic luminosity increasing downwards-higher observations are more sensitive to fainter galaxies. 
The only detections of possible DLA galaxies are shown by the two asterisks. The lowest curve is our 
measured sensitivity for Cycle 7; filled triangles represent direct imaging, while open triangles represent 
coronographic observations. The middle curve shows our predicted sensitivity in a 1300-second NICMOS 
integration, a half magnitude gain over our current limits. Finally, the upper curve is the predicted 
sensitivity of for "limiting" 10,000-second exposures, would push detection limits 2 mags fainter. 

There were only two plausible DLA detections and the best upper limits were H ~ 21.5 mag. When we 
account for the fact that more luminous galaxies are likely to have larger HI cross sections, we would have 
expected to detect 5 or 6 of the galaxies. If the DLAs were drawn from the general galaxy population, 
we should have detected 2 or 3 with 1.0L* or brighter, instead of zero. This discrepancy indicates 
that the luminosity function of DLAs is not consistent with the normal galaxy population. The objects 
associated with DLAs are on average fainter than normal galaxies we measure at either z=0 or 3. Deeper 
NICMOS/HST imaging can test the possibility that they are produced by dwarf galaxies such as the 
Magellanic Clouds. 

If deeper proposed NICMOS imaging of DLA systems continues to set even stricter upper limits on 
their continuum emission, it may be that DLAs trace a substantial baryon reservoir not associated with 
normal galaxies. 
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